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A New Measurement of the Stellar Mass Density at z~5: 
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We present a new measurement of the integrated stellar mass per comov- 

ing volume at redsliift 5 determined via spectral energy fitting drawn from a 

|7 j ■ sample of 214 photometrically-selected galaxies with z«^ nrp <26.5 in the so uth- 



(Ti ■ ern GOODS field. Following procedures introduced by lEyles et al.l (120051 ). we 

estimate stellar masses for various sub-samples for which reliable and uncon- 
CN ■ fused Spitzer IRAC detections are available. A spectroscopic sample of 14 of the 



most luminous sources with z = 4.92 provides a firm lower limit to the stellar 



10 . mass density of 1 x 10^ M© Mpc^^. Several galaxies in this sub-sample have 

t;;;-|- . masses of order 10^^ Mq implying significant earlier activity occurred in massive 

O , systems. We then consider a larger sample whose photometric redshifts in the 

O ■ publicly-available GOODS-MUSIC catalog lie in the range 4.4 < z < 5.6. Before 

^ . adopting the GOODS-MUSIC photometric redshifts, we check the accuracy of 



their photometry and explore the possibility of contamination by low-z galax- 
ies and low-mass stars. After excising probable stellar contaminants and using 
'^ . the ^^ssoLP " -^ color to exclude any remaining foreground red galaxies, we con- 

clude that 196 sources are likely to be at 2; ~ 5. The implied mass density from 



the unconfused IRAC fraction of this sample, scaled to the total available, is 
6 X 10^ Mq Mpc~^. We discuss the uncertainties as well as the likelihood that 
we have underestimated the true mass density. Including fainter and quiescent 
sources the total integrated density could be as high as 1 x 10^ Mq Mpc~^. Even 
accounting for 25% cosmic variance within a single GOODS field, such a high 
mass density only 1.2 Gyr after the Big Bang has interesting consequences for 
the implied past average star formation during the period when cosmic reion- 
ization is now thought to have taken place. Using the currently available (but 
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highly uncertain) rate of decline in the star formation history over 5 < 2; < 10, a 
better fit is obtained for the assembled mass at z ~ 5 if we admit significant dust 
extinction at early times or extend the luminosity function to very faint limits. 
An interesting consequence of the latter possibility is an abundant population of 
low luminosity sources just beyond the detection limits of current surveys. As 
mass density estimates improve at 2; ~ 5-6, our method is likely to provide one 
of the tightest constraints on the question of whether star forming sources were 
responsible for reionizing the Universe. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: starburst 
- galaxies: high redshift - ultraviolet: galaxies - surveys 



Introduction 



Finding the sources responsible for cosmic reionization is now the active frontier in stud- 
ies of galaxy formation. A number of independent arguments are focusing efforts on searches 
for star forming galaxies in the redshift interval 5 < z < 10. Studies of the optical depth in 
Lyman a absorption probed by high resolution spectra of the most distan t quasars suggest an 
upward transition in the neutral fraction beyond z ~ 5.5 (JFan et al.ll2006l ): these data suggest 
reionization was just ending at z ~ 6. In contrast, the optical depth of microwave photons 
to electron scattering derived from the angular power spect rum of the WMAP polarization- 
temperature cross-correlation function (ISpergel et al.ll2006l ) places a valuable upper bound 
on the reionization process corresponding to 2; ~ 10-20. 

Over the past several years, the quest to observe the most distant galaxies in the Universe 
has rapidly expanded to the point where the discovery of z ~ 5 — 6 star-forming galaxies 
has now become routine. Deep imaging surveys with the Hubble Space Telescope (HST) 
and 8-10 meter ground based telesco pes have u ncovered hundreds of galaxies at ^ :^ 5 
dlwata et al.lboosl: ISremer et all liooi) and z ~ 6 JSunker et al.lbooi IPickinson et al.ll2004J : 



Bouwens et al.ll2006l ) via the Lyman break galaxy (LBG) tec hnique pioneered by Steid el and 
collaborators to identify star-forming galaxies at z ~ 3 — 4 (ISteidel et al.lll996l . Il999l ). 



The consensus emerging from these studies, however, is that abundance of luminous 
galaxies is substantial l y less at 2: ^ 6 than a t z ~ 3 (IStanway et al.ll2003l : iBunker et al.l 12004 



Dickinson et al.ll2004l : iBouwens et al.ll2006l ). If this trend continues to fainter systems and 
higher redshift, then it may prove challenging to explain the earlier star formation activity 
necessary to fulfi ll reionization in th e redshift interval 5 < 2 < 10 implied by the quasar and 
WMAP studies (JBunker et al.ll2004l ). However, it has been suggested that the evolution in 
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the galaxy luminosity function between z = 3 and z = 6 is luminosity dependent: although 
the entire luminosity function is not yet well-constrained a.t z p^ 6, intrinsica lly fainter 
galaxies appear to become more abundant at earlier times (JBouwens et al.ll2006l ). If this is 
the case, then the bulk of reionizing photons could come from lower luminosity galaxies not 
yet adequately probed in deep surveys. 

As the redshift boundary of cosmic reionization narrows, so it becomes crucial to improve 
our understanding of the cosmic star formation history in the corresponding time interval. 
Unfortunately however, confirming even the most luminous sources in the range 7 < z < 10 
is challenging for current facilities. Although some candidate 2; ~ 7 — 10 galaxies have been 
identified in ACS and lensed surveys (Bouwens et al 2004b, Bouwens et al 2005, Richard 
et al 2006, Stark et al. 2007), these are generally too faint for spectroscopic study. The 
situation may not significantly improve for several years. 

This paper explores a more practical approach for constraining the amount of star 
formation prior to ^ ~ 5-6, namely the measurement of the integrated ste llar mass density 
at this epoch. Following the idea originally discussed by lStark &: Ellid (120061 ) . the stellar mass 
density at 2; ~ 5-6 must represents the integral of past activity. With adequate precision, such 
estimates can be used to independently verify the claimed decline in overall star formation 
to z ~ 10 and to assess whether the past activity is sufficient for cosmic reionization. 

The approach is made practical by the remarkable progress recently made in esti- 
mating stellar na asses at high redshift via the use of th e Infrared Ar r ay Ca mera (IRAC, 
Fazio et al.l |2004J ) onboard the Spitzer Space Telescope. lEgami et al.l (120051 ) first demon- 
strated the technique for one of the m ost distant k nown sources: a multiply-imaged pair 
with a photometric redshift of 2 ~ 6.8. lEyles et al.l (120051 ) later extended the technique for 
two spectroscopically-confirmed galaxies at z = 5.8, demonstrating the presence of massive 
galaxies (Msteiiar > 10^° Mq) with evolved stellar populations of ages ^ 100 Myr. 

The IRAC filters at 3.6 — 8.0 fim. probe the rest-frame optical at z ~ 5 — 6, providing a 
valuable indicator of established stellar populations and, indirectly, hinting at vigorous star 
formation activity at z > Q. Combining these data with deep broadband optical photometry 
from HST and 8-10 meter class ground based telescopes, spectral energy distributions (SEDs) 
can be compared with population synthesis models to constrain the age, star formation 
history and stellar ma sses of galax i es. T he initial discovery of massive (lO^^M©) galaxies 
at 2 ~ 6 p resented in Evles et al.l (120051 ) wa s subsequently c o nfirrn ed by the independent 
analysis of JYan et al.l (120051 ) . More recently, JMobasher et al.l (120051 ) identified a galaxy in 
the Hubble Ultra De ep Field (UDF) wi th a photometric redshift oi z ~ 6.5 (but see also 
the recent paper by iDunlop et al.l 120061 ) . If this high redshift is correct, then the MIPS 
and IRAC detections imply a very massive system of Msteiiar > 10^^ M0, providing further 
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evidence for significant star formation activity at z > 6 (jPanagia et al.ll2005l ) 



The studies of galaxy masses published thus far have focused on only a few individual 
systems. Although some studies (e.g. Stark & Ellis 2006) have attempted to infer the 
contribution of past star formation to cosmic reionization, without knowing how typical 
such massive galaxies are, it is difficult to make precise statements. Clearly what is needed 
is a census of the assembled stellar mass at high redshift. A comoving stellar mass density 
can be directly compared with various models of earlier star formation. 

In a companion paper, we compute the st ellar mass density at z ^ 6 from the z'-band 
dr opouts in GOOD S-South (JEyles et al.ll2006l ). A similar study of z'-drops was conducted 



m 



Yan et al.l (120061 ). However, the surface density of i'-hand dropout galaxies at z ^ 6 



with Spitzer detections is low. A more statistically-meaningful sample can be found using 
the z ^ 5 w-band dropouts. The age of the Universe at this time is only marginally older 
(1.2 Gyr c.f. 0.95 Gyr) yet larger, more representative, samples are available. In this paper 
we will examine the stellar mass density at z ~5 using sources to a limiting magnitude 
of ^8^0/ p ~ 26.5 selec ted from the Great Observatories Origins Deep Survey (GOODS, 
Giavalisco et al.ll2004bl ). We present an analysis of various subsamples at 2; ~ 5 drawn from 



a total of ^214 i;-band dropouts. The goal of the study is to establish whether the assembled 
stellar mass at 2; ^ 5 is consistent with current (and admittedly uncertain) estimates of the 
preceding star formation activity. If not, this might be taken to imply a significant component 
of star formation is missing, occurring either at lower intrinsic luminosities, obscured by dust, 
or at uncharted epochs {z > 10). 

A plan of the paper follows. In ^ ^ and §11 we introduce the various imaging and 
spectroscopic datasets, the photometric procedures and the selection of various subsamples 
of 2; ~ 5 galaxies. We describe the derivation of the stellar masses and comment on the 
uncertainties in §3 In §H1 we examine the implications for the star formation history at 
earlier times. 



W e adopt a cosmology consistent with the initial WMAP data release (jSpergel et al. 
2OO3I ) : a A-dominated, fiat universe with f^A = 0.7, ^m = 0. 3 and ifn = 70 /lynk ms"^ Mpc~^ 



All magnitudes in this paper are quoted in the AB system (lOke fc Gunnlll983l ) 



The GOODS-S Dataset 



In this paper, we continue our analyses of the Great Observatories Origins Deep Survey 
(GOODS). GOODS aims to bring together the most powerful space and ground-based fa- 
cilities to study the high-redshift universe across a wide range of wavelengths. We focus on 
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the southern GOODS field which has the greatest amount of multi-wavelength data essential 
for reliable stellar masses. The GOODS-S survey ar ea covers a total of 1 60 arcmin^ and is 
centered on the Chandra Deep Field South (CDF-S; iGiacconi et al.ll2002l ). 



2.1. A C5 Imaging 

Deep optic al imaging of GO ODS-S has been obtained with the Advanced Camera for 
Surveys f74C6'. iFord et al.l |2003| ) instrument onboard HST as part of a Treasury Program 
JGiavalisco et all l2nn4al )rThe Wide Field Camera on ACS has a field of 202x202 arcsec^ 
and a pixel scale of O'.'OS. The GOODS-South field was observed in the F435W (-B-band), 
F606W (w-band), F775W (SDSS-i') and F850LP (SDSS-^') broad-band filters for 3, 2.5, 2.5 
and 5 orbits, respectively over 16 pointings. 

Here we present an analysis of 2; ^ 5 galaxies making use of the publicly-available 
version-1.0 data-release of the ACS GOODS datall|. The reduced data have been 'drizzled' 
onto a large grid made up of 18 sections with a pixel scale of 0'.'03. Each section comprises 
an image of 8192 x 8192 pixels in size. 



2.2. Ground-Based Near-infrared Imaging 

Deep near-infrared observations of most of the GOODS-S field were obtained with the 
ISAAC camera on the Very Large Telescope (VLT) at the ESO Paranal Observatory as part 
of the ESO Large Programme: LP168.A-0485(A) (PI: C. Cesarsky). The publically available 
version- 1.5 data release includes 24 fully reduced ISAAC/ VLT pointings in the J and Kg- 
bandqj, covering ^ 160arcmin^. Additional details of the observations are to be presented 
in Vandame et al. {2006, in prep). The VLT images have a pixel scale of 0'.'15, a factor of 
five times larger than the drizzled ACS pixels. The median exposure times are 11.3 ksec in 
J, and 17.9 ksec in K^. 



^available from ftp : //archive . stsci . edu/pub/hlsp/goods/| 

^available from |h.ttp : //www . eso . org/science/goods/releases/20050930/ 
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2.3. Spectroscopy 



We also use publicly-available spectroscopy from the GOODS team to identify confirmed 
2; ~ 5 galaxies for futlier study. Multi-object spectroscopy was performed on the GOODS-S 
field with the F0RS2 instrument mounted at the Kueyen Unit Telescope of the VLT at 
ESO's Cerro Paranal Observatory as part of the ESO/G OODS Large Program LP1 70.A- 



0788 (PI Cesarsky). Details of the survey are presented in lVanzella et al.l (120021 . l2005l ). The 
primary selection criteria for placing objects on the slitmask was (^775^^ — ^85qlp) > 0-6 and 
^850LP < 25.0; objects with 0.45 < (^7751^ — -^^ssolp) < 0-6 were placed on the slitmask with 
lower priority. We make use of the VLT/F0RS2 spectroscopic catalogs from the version-2.0 
release which provide 725 unique redshift assignments with quality flags A, B, or C (where 
A=solid redshift, B=likely redshift, C=potential redshift). 



2.4. Spitzer Imaging 

Spitzer images of GOODS-S were obtained with the Infrared Array Camera (IRAC) and 
Multiband Imaging Photometry for Spitzer (MIPS) cameras on the Spitzer Space Telescope 
as part of the "Super Deep" Legacy programme (PID 169, Dickinson et al. in prep, Chary et 
al., in prep). The IRAC camera comprises four channels, each with a 256^ InSb array imaging 
a 5.2' X 5.2' field with a pixel size of ~ 1'.'22. Images were taken through four broad-band 
infrared filters, with central wavelengths at approximately Acent = 3.6 /xm, 4.5 /xm, 5.6 /xm 
and 8.0/im (channels 1-4), and widths of AAfwhm = 0.68, 0.87, 1.25, 2.53 yum respectively. 
The total exposure time in each channel is ~86ksec, depending on location. The data were 
taken in two epochs, with the telescope roll angle differing by 180°. In the first epoch, each 
filter covered a 10.0' x 10.0' area in GOODS-S; however, the area covered by channels 1 
and 3 (3.6 /zm and 5.6 /xm) was offset by 6.7 arcminutes from that covered by channels 2 
and 4 (4.5 /im and 8.0 /xm). Hence, only a portion of the GOODS-S field was observed in 
all 4 filters after the first epoch of observations. In the second epoch, the area covered by 
channels 1 and 3 in the first epoch was observed with channels 2 and 4 and vice versa. A 
central overlap region appeared in bo th epochs, and this deeper area intention ally contains 



the Hubble Ultra Deep Field (HUDF, ISeckwith et al.ll2003l : Bunker et al.l 120041 ) 



We analyze the publicly available Spitzer mosaics from the first and second epochs of the 
observations of GOODS-S q The data reduction pipehne employs a 'multidrizzle' technique 
similar to that used successfully on HST/ ACS GOODS data. This provides combined images 



•^available from^http: //data. spitzer. caltech.edu/popular/goods 
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with a pixel scale of 0'.'6. The magnitudes listed in this paper are determined from this 
'drizzled' data. We use the updated "Super Deep" epoch 1 images from the third data 
release (DR3) and the Super Deep epoch 2 images from the second data release (DR2). 



3. Photometric Samples 

The photometry we compute in this section will be used for two independent samples 
of 2; ~ 5 objects: a small sample of spectroscopically confirmed galaxies and a larger sample 
of photometrically selected galaxies. The spectroscopic sample will provide a robust lower 
limit to the 2; ~ 5 stellar mass density whereas the photometric sample will provide a more 
representative estimate of the integrated mass density. To obtain stellar masses of individual 
galaxies, we must have accurate photometry for both samples as well as photometric redshifts 
for the photometric sample. The reliability of the photometric redshifts is especially crucial 
since contamination by low-redshift interlopers could seriously skew our estimates of the 
total mass. 

We ob tain photometric red shifts from the GOODS MUSIC photometric catalog of 



GOODS-S (jGrazian et al.ll2006l ). This catalog uses uses 13-band SEDs from HST /ACS 
and Spitzer/YRKC photometry along with ground-based f/, J, & Ks to derive photometric 
redshifts. Before adopting the GOODS MUSIC photometric redshifts, we verify the accuracy 
of the photometry in the GOODS MUSIC catalog (discussed below) and test the reliability 
of their photometric redshifts which we discuss in §4.1. 

ACS photometry was obtained from the GOODS team rl.l catalog Q The photo- 
metric zeropoints adopted in the catalog on the AB magnitude system are 25.653, 26.493, 
25.641, and 24.843 for the i?435vy-band, Wgoew-band, i'tj^w band, and ZggQ^^p-band, respec- 
tively. We have corrected for the small amount of foreground Galactic extinction using 



the C05E/DIRBE & IRAS/lSSk dust maps of ISchleeel et al.l (119981 ): for the GOODS-S 
field, selective extinction is given by E(i? — V)= 0.008 mag. Magnitudes are measured in 
O'.'50-diameter apertures. Total magnitudes are derived from the aperture magnitudes by 
correcting for the small amount of light falling outside the aperture: 0.14 , 0.15, and 0.20 



magin the f606iy,'^775VK; ^^"^ ^ssoLP'^ands, respectively (jSirianni et al.ll2005l ). We note that 
GOODS website implies that the SExtractor parameter PHOT_APERTURES measures the 
radius of the photometric aperture, when it in fact measures the diameter. The correct 
interpretation has been applied to our dataset. 



available from http : //archive . stsci . edu/prepds/goods| 



Near-infrared photometry was performed with l"-diameter apertures using the ground- 
based near-infrared ISAAC images. The center of the photometric aperture was taken from 
the centroid of the GOODSrl.l catalog. The seeing varied across the ISAAC field as different 
tiles were taken over many nights, so we determined separate aperture corrections from 
unresolved sources for each tile. For the J- and Ks-hand images the seeing is typically good 
(FWHM = 0'.'4 — O'.'S), and the aperture corrections are ~ 0.3 — 0.5 mag, determined from 
bright but unsaturated isolated stars measured in 6"-diameter apertures. The 3a limiting 
AB-magnitudes in a l"-diameter aperture are J ~ 26.4 and Ks ~ 25.7, although these vary 
over the field because of different exposure times and seeing conditions. 

The details of the photometric analysis of t he Sp itzer images used in this paper are 



nearly identical to those presented in lEyles et al.l (120051 ). In order to maximize the signal-to- 



noise ratio (S/N) and minimize possible confusion with other foreground objects, we used 
a photometric aperture of diameter ^ 1.5 x FWHM for the IRAC images, appropriate for 
unr esolved objects (our compact souces are essentially unresolved at IRAC resolution, see 



e.g. iBremer et al.ll2004l ). The aperture diameters were 4, 4, 5 & 6 'drizzled' pixels for the 4 
channels (3.6, 4.5, 5.6 & 8.0 /xm), corresponding to 2'.'4, 2'.'4, 3'.'0, & 3'.'7. We used the IRAF 
digiphot .phot package to measure the enclosed flux at the coordinates determined by the 
ACS GOODS-rl.l catalogs, taking the residual background from an annulus between 12" 
and 24" radius. We applied aperture corrections to compensate for the flux falling outside 
the aperture: these were ~ 0.7 mag for the IRAC data, as determined from bright but 
unsaturated point sources in the images using large apertures. 

The noise for each of the four channels was checked in two different ways. First, we 
derived an estimate based on a Poisson model using the detector gain, number of frames com- 
bined, and the background counts (adding back the zodiacal background estimate subtracted 
by the pipehne but recorded in the header). Secondly, we measured the standard deviation 
in background counts of the images. As the mosaicking process introduces correlations be- 
tween pixels, we also made noise estimates using the individual pipehne basic calibrated data 
(BCD) images and assuming it decreased as the square root of the number of frames. These 
estimates lead to 3a limiting AB magnitudes of 26.5 and 26.1 using 2"4-diameter apertures 
in channels 1 and 2, respectively, and 23.8 and 23.5 in 3"0 and 3"7-diameter apertures in 
channels 3 and 4, respectively. There will be additional background fluctuations caused by 
faint galaxies (i.e. confusion noise), which will increase the noise. Both methods produce 
consistent estimates. 

The low spatial resolution of Spitzer results in frequent blending between nearby sources, 
making accurate photometry of individual objects difficult. We took great effort to ensure 
that objects in our sample were not contaminated by neighboring bright foreground sources. 
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We approach the IRAC contamination in shghtly different ways for the different subsamples 
of z~ 5 objects. Details are provided in §11 

We find that our photometry is consistent with that in the GOODS MUSIC cataog. 
The standard deviation between our photometry and the GOODS MUSIC photometry is 
0.13 mags, 0.03 mags, 0.13 mags for the z', J, and K band. This increases to 0.36 mags for 
the 3.6 micron IRAC photometry. 



4. Selection of z~5 Galaxies 



4.1. The Photometric Sample 

We make use of the e xtensive database of photometric redshifts in the publicly-available 
GOODS-MUSIC catalog (JGrazian et al.ll2006l ) to construct a sample of z ~ 5 candidates. De- 



tails o f the procedure used to compute the photometric redshifts are discussed in lGrazian et al 
(120061 ). As described below, we based our photometric selection on t he GOODS-MUSIC 



catalog rather than upon a more tradit ional v-band dropout technique (jBremer et al.l 12004 



Giavalisco et al.ll2004al : lYan et al.ll2005l ) on account of the improved performance in various 



tests. The principle difference is that the former method is based on fitting the entire SED. 

First, we consider the fidelity of the G OODS-MUSIC select i on of z ^ 5 galaxies with 
respect to the VLT spectroscopic results of IVanzella et al.l (J2002l . 120051 ). 21 galaxies within 
our 2 ~ 5 spectroscopic sample (see next section) have photometric redshifts in the GOODS- 
MUSIC catalog. 18 of these (>85%) have photometric redshifts in the 4.4 < z < 5.6 range 
with an average absolute scatter of < \zspec — Zphot\ >=0.07. Two of the three objects for 
which the photometric redshifts fail completely (e.g photometric redshifts of z ~ 1 — 2) have 
spectroscopic redshift quality grades of C: here it is possible that the photometric redshifts 
are actually correct. This test suggests the SED-fitting process is reasonably accurate. 

A further verification of the reliability of the photometric catalog concerns the implied 
rest-frame colors. Adopting a magnitu de limit of z'^^^rp < 2 6.5 (the 50% completeness 
limit for unresolved sources in GOODS iGiavalisco et al.ll2004al ). we find 214 objects with 
photometric redshifts between 4.4 < z < 5.6. Their rest-frame UV colors are in uniformly 
good agreement with those expected from the locus of star-forming galaxies at z ~ 5 (Figure 

We find that only 4 2% of the objects in ou r photometric catalog would have been 
selected in the traditional IGiavalisco et al.l (l2004al ) u-drop method. Examining the redshift 
tracks, it is clear that the v-drop method misses a significant fraction oi z = 4.5 — 5.5 
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star- forming galaxies (Figure [T]). This region of color-color space is not included in the 
traditional v-drop method to minimize the inclusion of low redshift contaminating galaxies. 
The GOODS-MUSIC photo-z sample (along with the criteria we impose below) represents 
an improvement to the traditional f-drop selection criteria as it takes the entire SED into 
consideration in assessing an object's redshift. 

While the GOODS-MUSIC photometric redshifts appear to be an excellent predictor 
of the true redshift, we remain vigilant to the possibility of a few catastrophic failures. The 
point is critical as the presence of any residual low redshift or stellar objects that are very 
bright at Spitzer wavelengths could lead to a significant overestimate of the stellar mass 
density. Recognizing there is a danger of removing true z ~5 sources, we conclude it is 
better to err on the conservative side. 

Because of their red colors, low- mass stars are a common contaminant of photometrically- 
located high redshift galaxy samples. Bright stars can be removed from high-z galaxy samples 
by selecting unresolved objects in the HST/ACS images. However, this technique begins to 
fail at fainter magnitudes as extragalactic objects may appear unresolved if observed at 
low S/N. Alternatively, stellar contaminants can be selected from our sample on the basis of 
their optical through near-infrared. We fit the SEDs of all objects in the photometric catalog 



with M, L, and T dwarf stellar templates (JLeggett et al.ll2002l : IWest et al.ll2005l : iKraus et al. 



20061 ). We construct a list of stellar contaminants by examining each object well-fit with 
stellar colors, only including sources without extended emission. Our final list consists of 11 
stars (5% of the total sample) with z'^^^j^p = 24.3 — 26.5, each of which we remove from our 
photo-z sample. 

Low-redshift galaxies with intrinsically red colors arising from dust extinction or an old 
stellar population commonly contaminate traditional dropout samples because their (vgoevi/- 
^775Vk) colors are similar to those of 2; ~ 5 star-forming objects. By considering the shape 
of the entire SED, low-z interlopers can often be identified and removed from high-redshift 
dropout samples. Since the GOODS-MUSIC photometric redshifts are computed using the 
entire SED, we expect the contamination rate from low-z galaxies to be low. Nevertheless, 
we believe it is important to explore the possibility that low-redshift galaxies may remain in 
the GOODS-MUSIC sample and examine the effects that possible contaminants may have 
on our final results. 

A simple way to estimate the contamination rate from low-z galaxies is to measure 
the rest-frame UV— optical colors of each of the objects in our sample. Unextincted star- 
forming objects at 2; ~ 5 typically have spectra that are roughly flat in f,^ (as a function 
of wavelength) between the Lyman break and rest-frame ~4000 A. In contrast, the colors 
of low-redshift contaminants are red in all filters. To quantify the expected difference in 
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rest-frame UV— o ptical colors between z ^ 5 sources and possible low-z contaminants, we 
examined a set of . 



Bruzual fc Charlotl (120031 ) population synthesis models. Elliptical galaxies 
at z ^1-2 with ages >2 Gyr have (^^gsoLp — <^) colors that vary between 1.4 - 1.6; whereas 
young (~ 100 Myr) star-forming galaxies at z ~5 with E(B— 1^)=0.0-0.2 have (zgsoLP ~ ^) 
colors ranging between —0.1 and —0.3 Accordingly, to test the low-z contamination rate, we 
adopt a (zgsoLP ~ J) > 1-0 threshold lI above which we consider galaxies to be possible low-z 
interlopers. Seven objects in our photo-z sample satisfy this color criterion. Six of the seven 
objects are relatively faint in the IRAC filters, and thus will hardly contribute to the total 
stellar mass of the sample. One of the objects (23_18055), however, is very bright in the near 
and mid-infrared (m3.e^m=21.1); if at z~5, its best-fit stellar mass would be 2x 10^^ Mq. 
Given that no objects are i dentified at z > 4 w ith stellar masses above 3x10^^ Mq in the 



0.8 sq deg UKIDSS survey (JDunlop et al.ll2006l ). we conclude that it is much more realistic 



to adopt a low-z interpretation for this object. So as not to bias our total mass estimates 
we remove the seven objects with (zggQ^p ~ J)> 1-0 from our photo-z sample, leaving 196 
objects. 

The final photometric sample is that for which the Spitzer IRAC images reveal a clear, 
unconfused, detection. Reliable stellar masses cannot otherwise be determined. We examined 
the Spitzer images of each of the 196 z ~5 candidates, classifying them as either (1) isolated 
and detected, (2) undetected, (3) confused or (4) hopelessly confused. In the subsequent 
analysis, we consider only those objects that are detected and isolated. Of the 196 candidates, 
72 are sufficiently uncontaminated to allow reliable estimates of the stellar mass. 

Table[T]lists the measured optical through infrared AB magnitudes (corrected to approx- 
imate total magnitudes through an aperture correction), colors, and photometric redshifts 
for the remaining 72 2 ~ 5 objects. 



4.2. The Spectroscopic Sample 

The F0RS2/VLT spectroscopic survey of the GOODS-S field identified 30 unique galax- 
ies in the 4.4 < z < 5.6 redshift range. The quality fiags associated with the redshift as- 
signments range from A (solid) to C (potential). As with the photometric sample, we adopt 
a magnitude limit of ^gsoLp < 26.5; this requirement excises one object (35_11820) from 
the sample. Given the possibility of uncertainties in the spectroscopic identification of those 
sources with C-grade redshifts, we examined their rest-frame ultraviolet colors (feoevF^^Wsvy) 
^s. (^775iy — -ZgsoLp) ^^ an additional criterion for selection (Figure [2]). 



^The precise value of this color discriminant is not critical in defining the final sample 
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Of the 29 remaining F0RS2 galaxies with spectrosco pic redshifts of 2; ~ 5, only 17 would 



be selected as f-drops using the iGiavalisco et al.l (l2004al ) selection criteria. An additional 8 



of the spectroscopically-confirmed z ~ 5 galaxies fall very near the f-drop selection window 
in the {vqqqw ~ '^775^) ^^- (^77511^ ~ ^ssolp) color-color plot. As their colors are consistent 
with the Bruzual & Chariot redshift tracks (plotted in Figure [2]) we include them in this 
sample. Three objects are apparently undetected in the ylCS" images of GOODS-S. Without 
the availability of rest-frame UV colors, we cannot confirm that the objects are truly located 
z ^5 via the presence of the Lyman break; we presume these were serendipitous detections 
and exclude them from the final spectroscopic sample. The final object (22_15184) is formally 
a B-dropout; its relatively blue (feoevy ~^775Vk) color is inconsistent with that expected from 
a v-drop. At the object's purported redshift {z = 5.08) Lya falls in the i'-band, making the 
intrinsic {vqqqw — i'rj^w) bluer than what is measured. Given the peculiar colors, we exclude 
it from the spectroscopic catalog. 

As before, we examined the Spitzer images of each of the 25 spectroscopically-confirmed 
galaxies for detections and the degree of confusion. These classifications are shown in Table [2l 
Five objects were isolated and detected in the Spitzer images, four sources were hopelessly 
confused, and the remaining 16 objects were marginally confused. For the 16 partially 
confused galaxies, we attempted to s ubtract the contr ibution from contaminating sources 



using the 'GALFIT' software package (jPeng et al.ll2002l ): this was deemed worthy given the 



need to maximize the information from the limited spectroscopic data. 

GALFIT constructs a two-dimensional model of the data according to specified input 
parameters (e.g., magnitude, position, axis ratio, effective radius), performs a convolution 
with the instrument point spread function (PSF), and fits the result to the data through 
an iterative x^ minimization process. We determined the PSF for each epoch and channel 
of the 'drizzled' Spitzer images by stacking 4 bright but isolated stars. For each galaxy we 
assumed a generalized Sersic surface brightness profile, where log/ oc r^/", and fit for the 
shape and index n. 

An automated script was developed to run GALFIT three times per source on a 12 x 
12 arcsec^ region surrounding the contaminated object for the IRAC images. In the first 
iteration, we held all source parameters fixed in the fitting process except the source mag- 
nitude, which was estim ated from the SExtractor source detection software version-2.2.1 



(JBertin fc Arnoutslll996l ) . All other input source parameters (e.g. position, axis ratio, posi- 
tion angle, effective radius, Sersic parameter) were estimated from a fit to the VLT i^^-band 
image. The higher spatial resolution of the Kg-hand allows better deblending and more 
accurate centroids to be derived for confusing objects in the IRAC images. In the second 
GALFIT iteration, we again determined input parameters using our fit to the i^^-band im- 
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age, but this time we allowed all parameters to vary. In the final iteration, we obtained 
the initial parameters by applying SExtractor to the IRAC channel 1 (3.6 yum) image and 
allowing all parameters to vary. For each source, we selected the most successful of the 
three GALFIT runs, on the basis of visual inspection of the residual image and the x^ value 
for the fit. Those sources (7 out of 16) for which GALFIT failed to satisfactorily subtract 
contaminating emission were removed from the sample (see Table [2]). The photometry of 
the remaining 14 galaxies are described in Table [31 



5. Stellar Mass Determination 

Although we have removed many sources from the original spectroscopic and photomet- 
ric samples, it is worth reminding that the degree of confusion in the IRAC images should, on 
average, be completely independent of the stellar mass of the z ~ 5 galaxy. Confusion in the 
IRAC images will normally arise from the overlapping isophotes of unrelated sources. Thus, 
if sources are believed to be at 2; ~ 5 on the basis of a spectroscopic redshift or the ACS and 
K-hand photometric SED, we can rescue a reasonable estimate of their contribution to the 
stellar mass density by scaling that determined for the unconfused sample using the relative 
numbers. 



5.1. Masses for the Spectroscopic Sample 

For those galaxies with confirmed spectroscopic redshifts, we estimate stellar masses 
by fitting population synthesis models to the observed SEDs. Applying this technique to 
galaxies for which the redshift is unknown may lead to significant uncertainty in the derived 
properties (JBundy et al.l l2005l : IShapley et al.ll2005l ). hence for the photometrically-selected 
sample, we infer stellar mass by applying the median mass-to-light ratio derived from the 
spectroscopic sample. 



We proceed as lEyles et al.l (120051 ) by fitting the latest iBruzual fc CharlotI (120031 ) stellar 
population s ynthesis models to the obs erved SEDs. We use the Padova evolutionary tracks 
preferred by iBruzual fc CharlotI (120031 ). The models utilise 221 age steps from 10^ to 2 x 
10^° yr, approximately logarithmically spaced. For each source, we d o not include ag e steps 
in excess of the age of the universe at z ^ 5 (~ 1.2 Gyr). Models with ISalpeterl (119551 ) initial 
mass fun c tions (IMF) were selected; although we also considered the effect of adopting a 
Chabrierl J2003h IMF. There are 6900 wavelength steps, with high resolution (FWHM 3 A) 
and 1 A pixels evenly-spaced over the wavelength range of 3300 A to 9500 A (unevenly spaced 
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outside this range). From the full range of metallicities offered by the code, we considered 
both solar and 0.25 Z© models. From several star formation histories available, a single 
stellar population (SSP - an instantaneous burst), a constant star formation rate (SFR), 
and exponentially decaying (r) SFR models with e-folding decay timescales r=10, 30, 50, 
70, 100, 200, 300, 500, 1000 Myr were used. 

For each of the galaxies in our sample, the filters were corrected to their rest-frame 
wavelengths by the appropriate redshift factor. The measured flux was folded through the 
filter transmission profiles, and the best-fit age model was computed by minimizing the 
reduced x^? using the measured errors on the magnitudes. The number of degrees of freedom 
is the number of independent data points (magnitudes in different wavebands) minus the 
number of parameters that we are fitting. The Bruzual & Chariot spectra are normalized to 
an initial mass of 1 Mq for the instantaneous burst (SSP) model, and an SFR of 1 Mq yr^^ 
for the continuous star formation model. The fitting routine returned the normalisation 
for the model which was the best-fit to the broad band photometry (i.e., minimized the 
reduced x^)- This normahzation was then used to calculate the corresponding best-fit total 
mass using the luminosity distance for the redshift of each source. When considering models 
other than an SSP (instantaneous burst), it was necessary to correct the total 'mass' values 
output by the fitting routine. For a constant SFR model, each of these masses needed to be 
multiplied by the corresponding best-fit age, since the B&C template normalization has the 
mass grow by 1 Mq yr~^. For the exponential decay models, the returned mass values were 
corrected by dividing by (1— e"*/"^), accounting for the decay timescale and the normalization 
of the B&C models (where M — ^ 1 Mq as t — > oo). The fits to the B&C models returned 
the 'total mass' which is the sum of the mass currently in stars, in stellar remnants, and 
in gas returned to the interstellar medium by evolved stars. For each best-fit model, we 
subsequently calculate the mass currently in stars for every galaxy, again using information 
from the B&C population synthesis code; we use this stellar mass in all future analysis. 

Although some of our data points (particularly from the HST/ ACS imaging) have 
S/N > 10, we set the minimum magnitude error to be A (mag) = 0.1 to account for calibra- 
tion uncertainties. Futhermore, we do not include data with photometric error above 0.72 
mags (1.5 cr). 

The presence of a strong spectral line in one of the broadband filters could significantly 
skew the SED fitting. Seven of the 15 galaxies in our spectroscopic sample show powerful Lya 
emission. Using the F0RS2 spectra, we compute and subtract the Lya contribution to the 
broadband flux; corrections rang e from 0.01 — 0.1 ma gs for most sources. Ha contamination 



could also be a significant issue - 1 Chary et al.l (120051 ) claim to find an excess due to Ha in a 



z=6.5 galaxy (in the 4.5/im band at that redshift). The sources in our sample are likely to 
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have Ha emission lines as well, which at z~5 fall in either the 3.6 /im or 4.5 yum IRAC filter. 
Indeed many of the SED fits discussed below show an excess at 3.6 yum. Without a direct 
measure of the Ha line strengths, we cannot robustly remove the line contamination. We 
estimate that, for most sources. Ha contributes ~ 10-20% of the measured broadband fiux 
by converting the inferred rest-frame UV star formation rate to an H-alpha luminosity via 
empirically-derived relations from Kennicutt (1998) assumin g < SFRunl SF Rjjy >=1.5-3 
due to dust extinction in agreement with observations, e.g. lErb et al.l (120031 ). To test the 
effects that Ha contamination may have on our sources, we re-fit all the objects in the 
spectroscopic sample, omitting the fiux information at 3.6 yum for objects with z < 5.2 and 
at 4.5 yum for objects with z > 5.2. We find that this does not significantly change the total 
stellar mass found in our spectroscopic sample. 

The de generacies associated with the derived best-fit parameters from SED fitting are 
well known (jShapley et al.ll200ll : iPapovich et al.ll200ll : IShapley et al.ll2005l ). The uncertain- 
ties primarily stem from a poor knowledge of the form of the star formation history, s ince 
the best fit age, dust extinction and star formation rate rely on this (jShapley et al.ll2005l ). In 
most cases, the data do not put strong constraints on the form of the star formation history; 
hence each fitted parameter typically has a range of values that produce acceptable fits. 

The inferred properties also rely on knowledge of the stellar initial mass function (IMF). 
There is little observational information constraining the IMF at high redshift. The spectrum 
of the z=2.7 gravitationally-lensed LEG cB58 appears to be inc onsistent with IMF s that have 
steep high-mass slopes or are truncate d at high stellar masse s (jPettini et al.ll2000l ). Whether 
this is typical among LBGs is unclear. iPapovich et al.l (120011 ) studied the effects that varying 
the IMF have on the best-fit parameters. Models with IMFs containing steep high-mass 
slopes (e.g. Scalo, Miller-Scalo) h ave redder integrat ed spectra and hence younger derived 
ages and lower extinction. As with lEyles et al.l (120051 ). we find good agreement between the 
properties inferred using a Salpeter IMF and a Chabrier IMF: the best-fit ages are nearly the 
same and the stellar masses are typically 30% lower when the Chabrier IMF is used. Here 
all masses are quoted for the Salpeter IMF in order to maintain consistency with previous 
estimates of stellar mass and star formation rates. 

In Figure HI we display the best-fit SEDs for each of the galaxies in our spectroscopic 
sample. The best-fitting model parameters are presented in Table HI The best fitting stellar 
masses of the galaxies range between 3x10^ Mq and 2x10^^ M©. Derived ages span three 
orders of magnitude, from 1 Myr to 1.1 Gyr, the age of the Universe at ^ ~ 5. Interestingly, 
three of our sources have stellar masses in excess of 10^^ Mq, values approaching the high 
stellar mass for the UDF source located by Mobasher et al (2005). Our results provides 
support for at least the presence of such galaxies even if their abundance remains uncertain. 
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Moreover, some of the less massive sources can only be fit wit h remarkably youn g ages (<20 



Myr) reminiscent of the lensed star forming source located by lEUis et al.l (120011 ) 



The total stellar mass of the subsample of spectroscopic galaxies is 5x 10^^ Mq. Clearly 
this estimate is an unrealistic lower limit to the total stellar mass since there are nearly twice 
as many objects known to be at 2; ~ 5 that certainly have stellar mass. 

Uncertainties on this mass arise from two main sources. First, the photometric error 
for each data-point in the SED translates into an uncertainty in the inferred stellar mass. 
Second, there is a range in acceptable masses that result from varying the age, extinction, 
and star formatio n history. With regard to the latter, we follow the approach outlined in 



Eyles et al.l (120051 ) where confidence intervals were explored for two sources. We present 
mass-age confidence intervals for two sources representative of our sample (Figure E]). Un- 
certainties in the inferred stellar mass of individual objects in Table 4 range between 30% 
and 50%. Objects detected at low signal to noise generally have larger uncertainties. Given 
the range of uncertainties, it seems reasonable to transfer a 50% uncertainty to all of our 
combined masses. 

The F0RS2 selection of galaxies was not geared specifically towards constructing a z~5 
sample; hence it is important to examine how the properties of spectroscopically-confirmed 
galaxies compares to the photometrically-selected sample. The median rest-frame UV color 
of the spectroscopic sample, (< ^/ysv^^-^ssoLP >)=0.17, is very similar to the photometrically- 
selected sample (< i^J^^r — ^ssqlp >=0.21). The key parameter for determining the stellar 
mass is the fiux in the IRAC filters. The histogram of IRAC 3.6 /im fiuxes for the spectro- 
scopic and photometric samples is given in Figure [31 The spectroscopic sample does contain 
a larger fraction of Spitzer-bright (e.g mae/^m < 23) objects, but this is reasonable if the 
overall rest-frame mass/light distribution is fairly similar across the population. 



5.2. Masses for the Photometric Sample 

To estimate stellar masses for the photometric sample we compute the best-fitting rest- 
frame l^-band mass to light ratio of each galaxy that is unconfused in the IRAC images and 
multiply by the luminosity derived from the IRAC fiux. 

The best-fitting M/Ly is determined for each galaxy from its ^^ssolp ~ ''^ssfim color 
(corresponding to the ratio of rest-frame UV and optical fiuxes). If we assume the typical 
galaxy in rest-frame UV selected samples at z^5— 6 has little dust as seems reasonable (see 
Table 4 and Eyles et al. 2006), then the ^ssolp ~ ^^^'3.6^lm color is correlated with the age of 
the galaxy and hence, for a given IMF and star formation history, its M/Ly ratio. 
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This is done for a given galaxy in the sample by first computing -ZgsoLp— ""^3.6^™ colors for 
Bruzual- Chariot models (redshifted to the galaxy's photometric redshift) with ages ranging 
between and 1.2 Gyr (roughly the lookback time at the redshift of the galaxy). We then find 
the model that produces the ^^ssolp "^s-e/xm color closest to that observed for a given galaxy. 
This model is taken to have the "best-fit" age and M/Ly ratio for this particular galaxy. At 
2; ~ 5, the 3.6/im IRAC filter covers the rest-frame V-band; hence, we convert the 3.6 /im fiux 
to a luminosity (assuming z = Zphot) and multiply it by the best-fitting M/Ly to compute 
the stellar mass. For each galaxy, best-fitting stellar masses are computed for the same range 
of single-component star formation histories used to fit the spectroscopic sample. The stellar 
mass we assign to each galaxy is taken from the star formation history that produces the 
best-fitting z'^^^i^p — m^,Q^j_m colors. We obtain an estimate of the systematic uncertainty in 
the mass by considering the range inferred from the different star formation histories and 
ages that provide a good fit (e.g. Ax^ = X^ ^ Xmin < 1)^*^ ^^^ observed -ZgsoLP ^ ""^s.e^m 
color. 



We note that the observed 3.6/im luminosity is not equivalent to a rest-frame V-band 
luminosity for all redshifts. The 3.6/im band shifts between rest-frame 5500-6700 A for 
z = 4.4-5.5. To test the systematic offsets introduced by relying on M/Ly to derive masses, 
we compare the mass of the spectroscopic sample derived in the manner described above to 
the mass from SED fitting. We find the median offset between the two methods is 40%. 

The total stellar mass extracted from the 72 z ~ 5 sources that are uncontaminated 
in the Spitzer images is 5 — 9x10^^ Mq with a best-fit value of 7x 10^^ Mq. The median 
stellar mass in the sample is 6x10^ M©. If we make the reasonable assumption that the 
distribution of stellar masses is independent of IRAC contamination, we can estimate the 
stellar mass in IRAC contaminated galaxies by multiplying the stellar mass derived from 
the uncontaminated z ^ 5 galaxies by the ratio of the total number of 2; ~ 5 sources to the 
number of uncontaminated z ^ 5 sources. Following this reasoning, the total stellar mass 
for the photometric sample becomes 2xlO^^M0. Taking the full range of single-component 
star formation histories into consideration, this total stellar mass could lie between 2 and 
3x 10^2 Mq. 



5.3. Comoving Mass Densities 

To derive the comoving stellar mass densities from the above totals, we need to estimate 
the redshift-dependent selection function in the 160 arcmin^ GOODS-S field between 4.4 < 
z< 5.6. Although the total possible comoving volume is 5.6 xlO^ Mpc^, the effective volume 
is less than this value due to sample incompleteness arising as a result of objects being 
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scattered faintward of the magnitude limit or out of the color-selection window. 

In order to account for these lumin osity and red s hift b iases, we compute an effective 
survey volume following the approach of ISteidel et al.l (119991 ) using 



Ves{m) = / dzp{m,z) — 



where p{m, z) is the probabihty of detecting a galaxy at redshift z and apparent z' magnitude 
m, and dz ^ is the comoving volume per unit solid angle in a slice dz at redshift z = 4.4—5.6. 

We compute the probability function p{m, z) by putting thousands of fake galaxies into 
the GOODS images and recreating a photometric catalog for the new image using identical 
selection parameters used in generating the GOODSvl.l catalogs. The apparent magnitudes 
of the fake galaxies span -ZgsoLp = 22 — 27 in steps of A?7i=0.5 and redshifts ranging between 
z = A and 2; = 6 in steps of Az=Q.l. The sizes of the fake galaxies are consistent with 
distriution of half-light radii derived for 2; ~ 5 galaxies in Ferguson et al. (2004). The 
colors of the fake galaxies depend on the galaxy redshift and SED. We adopt the SED of a 
Bruzual-Charlot model with constant star formation history, an age of 100 Myr, and no dust 
as the intrinsic rest-frame SED of the fake galaxies. Allowing for a selective extinction of E(B- 
V)=0.1 in the fake galaxies' SEDs decreases the effective volume by roughly 5%, which would 
not significantly change our final mass density estimates. The colors are then determined 
for galaxies at each redshift in a manner similar to that which we described in §5.1. The 
probability function, p{m, z) is then given by the fraction of fake galaxies with apparent 
magnitude, m, and redshift, z, that are brighter than the magnitude limit and satisfy the 
dropout color selection criteria. Since our selection is based on photometric redshifts, we 
adopt color-criteria that are appropriate for our photometric sample (feoew/ ~ "^7751^ > 0-9 

and «775p^ — ^850LP > i--o). 

The effective volume probed is 5.2x10^ Mpc^ at z'^^qlp = 23 and z = b where we are 
nearly 100% complete and falls to 1.2x10^ Mpc^ at Zg5Qj^p=26.5. The stellar mass density 
inferred from our z ~5 candidates is thus 5 — 8x10^ Mq Mpc^^, with a best-fit value of 6x10^ 
Mq Mpc~^. The robust lower limit from our spectroscopic sample is 1x10^ Mq Mpc~^. 

Our inferred stellar mass density is most likely an underestimate of the total value at 
z ~ 5 for several reasons. Foremost, the survey is only sensitive to the most luminous and 
perhaps most massive galaxies since we only included objects with significance above 3cr at 
3.6/im. 

Second, an additional reservoir of stellar mass may be contained in objects that are 
not currently forming stars and hence are very faint in the rest-frame ultraviolet. At 
z ~ 3, LBGs contribute only 17% of the stellar mass density in the most massive sources 
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( Ivan Dokkum et al.ll2006l ): the remaining fraction is likely contained in objects that are not 
actively forming stars. Although this fraction of quiescent sources is probably much lower 
at earlier times, we conservatively estimate the total mass density could rise by a further 
factor of two. 

In summary, therefore, we derive a firm lower limit to the stellar mass density at 2 ~ 
5 of 1x10^ Mq Mpc~^, a reasonable estimate of the total observed population of 5-8x10^ 
Mq Mpc""^ and cannot exclude undetected sources which would increase the total to 1x10^ 
M0 Mpc~'^. Although the overall estimates span a factor of 2, we emphasize that the 
spectroscopic sample is clearly a significant underestimate of the observed population. 



6. Implications for the Previous Star Formation History 

In the foregoing we have attempted to put the first bounds on the stellar mass density 
at 2 ~ 5, 1.2 Gyr after the Big Bang and about 800 Myr after z ~ 10. We emphasize that 
there are considerable uncertainties in the various steps in our analysis. First, to derive 
stellar mass estimates, we had to cull our samples to those with reliable IRAC detections, 
later scaling on the assumption that they represent a fair subset of the spectroscopic and 
photometric populations. For our spectroscopic sample, our fitting procedure gives mass 
estimates that span a wide range depending on the assumed star formation history. Finally, 
we assumed a median visual mass/light ratio for the photometric sample derived from that 
for the spectroscopic sample. 

Probably the dominant error in deriving the total mass density is not the numerical scal- 
ing factors , but rather the intrinsic uncertainty in estimating the masses of individual galax- 
ies. D etailed work at lower redshift (IShapley et al.ll2005l : iBundy et al.ll2005l : iPapovich et al. 
2OO5I ) has shown that inferences of the stellar mass from SED fitting yield mass estimates 
that contain typical uncertainties of ~30%. The errors certainly increase slightly when con- 
sidering objects at higher redshift; however, our error estimates (§5.1) suggest that the stellar 
mass estimates of objects in our spectroscopic sample are typically 50%, possibly more. 

Notwithstanding the uncertainties, it is interesting to now consider the implications of 
our derived mass density. The star formation rate density (SFRD) of bright (> 0.3 L*^3) 
star- forming galaxie s at z ^ 5 — 10 appears to d ecline continuously toward higher redshift 



(JBunker et al.l l2004l : iBouwens et al.l l2004j . |2005| ). However, current observations may be 
missing a substantial fraction of star formation either because it is enshrouded in dust, 
too faint to be detected with current facilities, or located at redshifts uncharted by current 
telescopes {z > 10). Comparing the co moving density of assembled stellar mass at z c^ 5 
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with estimates derived from models of the previous star formation history enables us to test 
these possiblities, thereby providing constraints beyond direct reach of current facilities. 



Taking data on the SFRD from the recent li terature (JGiavalisco et al.ll2004al : iBouwens et al. 



2004 : IStanwayll2004l : iBouwens et al.ll2005l . l2006l ). we fit the redshift dependence with a simple 
functional form over 2; ~ 5 — 10. In cases where the SFRD was not evaluated down to the 
adopted fiducial luminosity limit (O.IL* at z=3), we compute the additional contribution by 
integrating the luminosity function assuming the Schechter function parameters derived in 
each paper. Since the data at z > 7 do not allow for the robust derivation of the form of 
the luminosity function, we assume the shape of the luminosity function remains constant 
before 2; ~ 6. 



As lBouwens et al.l ( 120061 ) discuss, at z ~6 there is some disagreement in the value of the 
SFRD. The disagreement stems primarily from whether the shape of the luminosity function 
is evolving. Bouwens et al argue for a decrease in the characteristic luminosity and an increase 
in the faint end slope prior to z~3, so we adopt the redshift-dependent LF parameters derived 
in Bouwens et al. (2006) and Stanway (2005) and integrate accordingly. We find that the 
SFRD (integrated down to 0.1 L^^g) can be fit reasonably well by p oc (1 + z)~^-^ between 
z=5-10 (Figure ED. 

The stellar mass density obtained by integrating this function over time between z~10 
and 2; ~ 5 is lower than that derived from the photometric sample of 2; ~ 5 objects in this 
paper (Figure [7]). We note that the integral of the star formation rate density as a function 
of redshift overestimates the mass density. This is because we do not account for the mass 
that is returned to the interstellar medium in stellar winds and stellar deaths. This can 



be quantified using the iBruzual fc CharlotI (120031 ) software (we do this to compute stellar 
masses, see §5.1), but it is a complicated function of the average star formation history and 
age, which are not well-constrained. We find that this effect could reduce the stellar masses 
inferred from integrating the observed star formation rate densities by up to ~30% (assuming 
a 1 Gyr instantaneous burst). Given the uncertainties, we do not adjust the curves in Figure 
7 by this factor, but we note that this effect further enhances the discrepancy between the 
observed stellar mass density at z^ 5 and that which can be accounted for by previous star 
formation. Therefore, the observed stellar mass of the 2; ~ 5 galaxies in GOODS-S either 
implies a significant amount of dust extinction or that not all star formation at z>5 has been 
observed in current surveys. 

To examine the amount of star-formation that may be hidden in low-luminosity sys- 
tems, we integrate the luminosity fu nctions to zero lumin osity utilizing faint-end slopes of 
a. = —1.73 (as measure d at z~6 in iBouwens et al.l l2006l ) and a = —1.9 (as suggested by 
Yan Sz WindhorstI |2004| ) and integrate the luminosity function to zero luminosity. 



-21- 



Star formation is unlikely to occur at very low luminosities because of radiative feedback 
processes and (after reionization) a photoionized IGM which raises the cosmological Jeans 
mass. We nonetheless take this extreme approach to place an upper limit on the amount of 
unextincted star formation. Assuming no evolution in the shape of the LF between z ^ 5 
and 2; ~ 11, this increases the predicted stellar mass at z ~ 5 by an extra factor of 2.3 for 
faint end slopes of a = —1.73; if we instead consider an extreme faint-end slope of a = —1.9, 
the star formation rate density increases by a factor of 6.8. This gives a better account of 
the assembled mass and if correct has i nteres ting consequences for higher redshift surveys 
probing to low luminosities IStark et al.l (120071 ). 



A significant amount of star formation may also be enshrouded by dust. However, 
recent observations have shown that the rest-frame UV slope of z ~6 galaxies is actually 
somewh at bluer than that at z :^3 suggesting t hat the mean dust ex tinction declines between 
z ~3-6 (IStanway et al.ll2005l : lYan et al.ll2005l : iBouwens et al.ll2006l ). Taking the empirically 
derived fit relating the extinction at 1600 A (Aieoo) to the UV slope /5, Ai6oo=4.43 + 1.99/? 
(JMeurer et al.lll999l ) yields an overall attenuation factor of x ~1-1.5 at z ~6 for two different 



estimates of t he UV continuum slo pe at z 
f3 = -2.0 from lBouwens et al.l[2006[ ). 



6 (/? = —2.2 from IStanway et al.l l2005l and 



Hence, the expected extinction correction to the SFR density at z ~ 5 — 10 could in 
principle account for the stellar mass contained in the photometric sample in this paper if the 
Bouwens et al.l (120061 ) estimate of the UV continuum slope is correct. If there exists either a 
significant population of quiescent massive galaxies or low-mass galaxies below the 3a 3.6 fun 
flux limit imposed on the data, a significant amount of low-luminosity star-forming galaxies 
would be required to assemble the stellar mass. Future studies will test this hypothesis. 



Conclusions 



We have argued that the assembled stellar mass density at high redshift provides a 
valuable constraint on the past star formation history and, with improved precision, may 
ultimately indicate whether there was sufficient star formation in the previous ~500-900 Myr 
to reionize the intergalactic medium. 

We have demonstrated both the promise and limitations of this method by computing 
the co moving stellar mass density at z ~ 5. Following the ideas discussed in IStark fc Ellis 
(I2OO6I ) , we use the stellar mass density to constrain the amount of star formation at earlier 
times. 



We detail our findings below. 
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1. We construct a sample of 25 spectroscopically confirmed z ~5 objects in GOODS- 
S (14 of which are uncontaminated in the Spitzer data) to place a robust lower limit on 
the comoving stellar mas s density. Fitting the SEDs of these objects to templates from 



Bruzual fc CharlotI (120031 ) populations synthesis models, we infer a total comoving stellar 



mass density of 1x10^ M0 Mpc ''. 

2. We construct a sample of 2 ~ 5 galaxies using the photometric redshifts of the 
GOODS-MUSIC catalog. After removing likely stellar and low-z contaminants, 196 objects 
remain in the sample. Computing the stellar mass from the 72 objects that are uncontam- 
inated by nearby sources in the Spitzer data, we estimate a stellar mass density of 6x10^ 
Mq Mpc~^. Systematic uncertainty in the star formation history causes this value to be 
uncertain at the 30% level. 

3. The total comoving stellar mass density (6x10^ Mq Mpc^^) represents a lower limit 
for several reasons. First, robust stellar mass estimates are only attainable for reasonably 
massive galaxies; hence the estimates presented in this paper do not include the contribution 
from low-mass systems. Second, we require objects to be bright in the rest-frame UV (and 
hence actively forming stars) for selection into our sample. If there is a large population of 
quiescent galaxies at z ~ 5, the total stellar mass density may be significantly higher than 
estimated. We estimate that the stellar mass density of massive galaxies is unlikely to exceed 
1x10^ Mq Mpc-3. 

4. The estimated comoving stellar mass density at 2; ^ 5 suggests that current obser- 
vations may be missing some star formation at z > 5. The missing star formation could, 
however, be accomodated by extincted star formation in LBGs currently seen at 2; ~ 6-10 or 
in low-luminosity star-forming systems below the detection threshold of current observations. 
In the latter case, our results have important implications for searches for low luminsoity SF 
systems at high redshift. 
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Table 1. Photometric catalog of 2; ~ 5 galaxies in GOODS-S Field 



ID 



RA (,T2000) Doc (J2000) 



Ks 



3.6 /^m 



4.5 /im 



Zphot 



44.2919 


03 32 9.054 


-27 43 51.85 


27.58±0.16 


26.02±0.08 


25.65±0.07 


25.37±0.26 


25.50±0.44 


24.01±0.09 


24.29±0.14 


4.590 


42.3601 


03 32 


10.64 


-27 50 29.15 


26.83±0.08 


25.69±0.06 


25.60±0.07 


26.04±0.42 


25.37±0.37 


24.93±0.20 


25.41±0.43 


4.480 


33.4001 


03 32 


11.44 


-27 47 38.63 


27.88±0.18 


25.95±0.06 


25.77±0.07 


25.43±0.31 


25.46±0.71 


24.00±0.09 


24.31±0.17 


4.720 


33.4496 


03 32 


12.42 


-27 47 2.483 


27.83±0.20 


26.35±0.10 


26.21±0.12 


26.35±0.59 


26.03±0.98 


>25.8 


>25.4 


4.480 


33.4687 


03 32 


12.78 


-27 48 2.599 


27.47±0.11 


26.15±0.07 


25.96±0.08 


25.48±0.26 


25.58±0.62 


24.29±0.10 


24.95±0.23 


4.590 


34.4915 


03 32 


13.25 


-27 43 8.289 


27.90±0.22 


26.32±0.10 


26.16±0.11 


25.93±0.41 


25.85±0.68 


25.59±0.41 


>25.4 


4.680 


35.5207 


03 32 


13.88 


-27 41 48.54 


29.35±0.83 


27.27±0.23 


26.49±0.15 


29.30±9.44 


>25.3 


25.88±0.61 


>25.4 


5.340 


33.5533 


03 32 


14.49 


-27 49 32.69 


26.67±0.07 


25.69±0.06 


25.53±0.06 


26.21±0.56 


25.53±0.48 


25.07±0.33 


25.85±0.72 


4.460 


33.5986 


03 32 


15.35 


-27 49 36.08 


27.78±0.19 


26.15±0.09 


26.35±0.13 


25.79±0.38 


25.40±0.41 


25.27±0.58 


26.18±0.44 


4.610 


33.6438 


03 32 


16.17 


-27 46 41.59 


28.74±0.47 


26.25±0.10 


26.00±0.10 


26.63±0.75 


>25.3 


25.16±0.22 


26.04±0.64 


5.080 


33.6440 


03 32 


16.17 


-27 48 19.42 


27.83±0.20 


26.39±0.11 


26.25±0.13 


25.71±0.33 


25.00±0.38 


24.88±0.15 


25.28±0.31 


4.700 


33.6519 


03 32 


16.34 


-27 48 31.99 


27.87±0.21 


26.06±0.08 


26.13±0.11 


>26.0 


>25.3 


25.79±0.31 


26.19±0.68 


4.490 


33.6575 


03 32 


16.45 


-27 46 39.24 


29.45±0.89 


26.49±0.12 


26.11±0.11 


25.90±0.39 


25.08±0.41 


24.22±0.09 


25.00±0.25 


5.210 


32.6854 


03 32 


16.98 


-27 51 23.17 


27.41±0.14 


25.62±0.05 


25.70±0.08 


25.97±0.44 


24.48±0.19 


23.97±0.07 


24.44±0.16 


4.550 


35.6867 


03 32 


17.00 


-27 41 13.71 


26.89±0.08 


25.38±0.04 


25.13±0.04 


24.92±0.23 


24.37±0.24 


23.43±0.04 


23.82±0.08 


4.590 


32.8020 


03 32 


18.91 


-27 53 2.746 


27.77±0.19 


25.13±0.03 


24.49±0.03 


24.74±0.13 


24.06±0.13 


22.73±0.02 


22.74±0.03 


5.550 


31.8593 


03 32 


19.96 


-27 54 58.98 


28.64±0.56 


26.84±0.21 


25.83±0.11 


25.50±0.34 


25.59±0.45 


24.66±0.11 


25.23±0.26 


5.320 


31.9014 


03 32 20.70 


-27 55 36.14 


26.70±0.10 


25.67±0.08 


25.72±0.10 


25.99±0.54 


24.99±0.26 


24.08±0.07 


24.32±0.12 


4.520 


33.9184 


03 32 21.01 


-27 49 59.16 


27.23±0.12 


26.03±0.08 


26.16±0.11 


26.67±0.63 


>25.3 


>25.8 


>25.4 


4.580 


33.9338 


03 32 21.28 


-27 49 59.67 


28.54±0.38 


26.84±0.16 


26.14±0.11 


25.85±0.30 


26.64±0.92 


25.76±0.26 


25.35±0.26 


5.500 


33.9677 


03 32 21.82 


-27 50 3.346 


28.74±0.46 


26.49±0.12 


26.48±0.15 


>26.0 


25.48±0.31 


25.56±0.28 


>25.4 


4.800 


34.9738 


03 32 21.93 


-27 45 33.07 


28.22±0.29 


26.20±0.09 


25.82±0.09 


26.20±0.51 


24.90±0.23 


24.31±0.09 


24.72±0.18 


4.800 


33.9812 


03 32 22.02 


-27 46 42.89 


26.76±0.08 


25.41±0.05 


25.26±0.05 


25.94±0.48 


24.84±0.25 


23.84±0.06 


24.05±0.10 


4.510 


34.9822 


03 32 22.03 


-27 45 29.31 


27.38±0.14 


26.41±0.11 


26.28±0.13 


26.00±0.42 


>25.3 


25.30±0.22 


25.42±0.33 


4.570 


33.10064 


03 32 22.44 


-27 47 46.17 


28.58±0.40 


26.64±0.14 


26.34±0.13 


26.46±0.77 


25.57±0.47 


24.86±0.13 


25.22±0.27 


5.020 


32.10232 


03 32 22.71 


-27 51 54.40 


27.90±0.25 


26.14±0.10 


25.68±0.08 


25.58±0.28 


25.03±0.24 


24.27±0.08 


24.82±0.16 


5.050 


33.10340 


03 32 22.88 


-27 47 27.56 


26.64±0.07 


24.94±0.04 


24.84±0.04 


24.55±0.13 


24.59±0.16 


23.75±0.05 


24.01±0.10 


4.440 


31.10974 


03 32 24.00 


-27 54 59.79 


27.51±0.16 


25.46±0.05 


24.73±0.03 


24.56±0.16 


24.89±0.25 


25.39±0.41 


25.88±0.66 


5.380 


32.11635 


03 32 25.02 


-27 50 24.49 


29.10±0.56 


27.13±0.18 


26.05±0.09 


25.61±0.24 


25.42±0.30 


24.68±0.13 


24.67±0.17 


5.430 


33.13701 


03 32 27.94 


-27 46 18.57 


26.37±0.06 


25.18±0.04 


25.22±0.05 


25.09±0.21 


24.44±0.15 


24.02±0.10 


24.13±0.13 


4.480 


34.14195 


03 32 28.70 


-27 42 28.95 


28.11±0.21 


26.23±0.08 


26.03±0.08 


26.82±0.84 


>25.3 


25.67±0.34 


>25.4 


4.840 


23.15316 


03 32 30.28 


-27 49 22.01 


28.04±0.25 


26.47±0.12 


25.88±0.09 


26.30±0.47 


>25.3 


>25.8 


>25.4 


5.200 


22.15851 


03 32 31.07 


-27 51 17.85 


28.91±0.54 


26.17±0.09 


26.04±0.10 


26.06±0.40 


24.89±0.22 


24.62±0.17 


24.93±0.25 


4.820 


23.16055 


03 32 31.37 


-27 48 13.81 


28.08±0.26 


26.20±0.10 


25.89±0.09 


25.09±0.20 


25.66±0.39 


24.94±0.14 


25.42±0.32 


4.990 


22.17535 


03 32 33.69 


-27 53 21.62 


29.06±0.61 


27.28±0.24 


26.39±0.14 


>26.0 


>25.3 


26.28±0.60 


>25.4 


5.360 



to 
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ID 


RA (.72000) 


Doc (J2000) 


V 


i 


z' 


J 


Ks 


3.6 fim 


4.5 ^m 


Zphot 


23.17728 


03 32 33.98 


-27 48 2.043 


27.61±0.17 


26.10±0.09 


25.90±0.09 


25.74±0.29 


24.52±0.13 


24.08±0.07 


24.24±0.11 


4.470 


23.18716 


03 32 35.45 


-27 49 35.20 


29.41±0.60 


26.63±0.10 


26.35±0.09 


25.71±0.28 


25.94±0.49 


25.12±0.22 


25.91±0.56 


4.930 


22.19011 


03 32 35.89 


-27 52 44.02 


27.78±0.19 


26.51±0.12 


26.27±0.12 


26.92±0.85 


>25.3 


25.70±0.34 


>25.4 


4.740 


24.19118 


03 32 36.08 


-27 44 3.942 


27.28±0.12 


26.07±0.08 


25.75±0.08 


24.80±0.12 


25.22±0.26 


24.38±0.13 


25.01±0.27 


4.550 


23.19268 


03 32 36.30 


-27 49 52.79 


27.44±0.14 


26.03±0.08 


25.86±0.09 


25.36±0.21 


24.70±0.16 


23.96±0.09 


24.56±0.18 


4.510 


24.19435 


03 32 36.49 


-27 43 53.46 


28.52±0.37 


26.93±0.17 


26.43±0.14 


25.83±0.31 


25.93±0.55 


25.01±0.19 


25.59±0.42 


4.620 


25.19912 


03 32 37.25 


-27 42 2.570 


28.47±0.35 


26.42±0.11 


26.00±0.10 


25.81±0.32 


>25.3 


25.04±0.18 


25.53±0.41 


5.180 


22.20159 


03 32 37.62 


-27 50 22.38 


>29.5 


27.18±0.22 


26.24±0.12 


25.37±0.20 


25.54±0.36 


24.64±0.14 


24.66±0.17 


5.510 


22.20304 


03 32 37.86 


-27 52 59.10 


27.43±0.15 


26.47±0.12 


26.25±0.13 


>26.0 


>25.3 


25.68±0.45 


25.92±0.57 


4.520 


23.20360 


03 32 37.95 


-27 47 11.05 


27.39±0.13 


25.80±0.07 


26.21±0.12 


>26.0 


>25.3 


24.93±0.19 


>25.4 


4.690 


22.21669 


03 32 40.08 


-27 50 49.60 


27.72±0.18 


26.28±0.10 


26.05±0.10 


25.74±0.30 


25.56±0.40 


25.45±0.29 


26.17±0.64 


4.520 


24.22091 


03 32 40.85 


-27 45 46.25 


28.03±0.18 


26.14±0.07 


25.33±0.04 


25.12±0.14 


25.22±0.32 


24.76±0.12 


24.45±0.14 


5.400 


23.22354 


03 32 41.34 


-27 48 43.13 


28.02±0.24 


26.31±0.11 


26.14±0.11 


>26.0 


26.02±0.62 


24.72±0.13 


25.46±0.33 


4.530 


25.22925 


03 32 42.36 


-27 41 14.87 


27.38±0.21 


25.94±0.13 


25.41±0.09 


26.10±0.62 


24.65±0.36 


24.12±0.09 


24.77±0.21 


5.010 


24.23215 


03 32 42.95 


-27 43 39.65 


29.16±0.67 


26.68±0.14 


26.09±0.11 


>26.0 


>25.3 


25.94±0.59 


25.84±0.63 


5.200 


24.23395 


03 32 43.30 


-27 43 10.59 


28.73±0.45 


26.14±0.09 


26.25±0.12 


25.60±0.32 


>25.3 


24.74±0.16 


25.10±0.30 


4.780 


23.23515 


03 32 43.53 


-27 49 19.21 


27.89±0.21 


25.80±0.07 


25.53±0.07 


26.05±0.34 


25.84±0.50 


24.89±0.15 


25.92±0.49 


5.090 


22.25323 


03 32 47.58 


-27 52 28.18 


28.49±0.34 


26.43±0.11 


26.24±0.11 


25.43±0.26 


>25.3 


25.43±0.29 


>25.4 


4.760 


13.25544 


03 32 48.14 


-27 48 17.69 


27.95±0.23 


26.04±0.09 


25.41±0.06 


25.64±0.28 


24.89±0.25 


24.65±0.12 


24.97±0.24 


5.170 


14.25620 


03 32 48.33 


-27 45 38.90 


27.86±0.21 


26.24±0.11 


26.37±0.14 


26.42±0.49 


25.75±0.52 


>25.8 


>25.4 


4.610 


12.25696 


03 32 48.53 


-27 54 25.67 


27.44±0.12 


26.28±0.08 


26.24±0.10 


26.49±0.77 


>25.3 


24.89±0.20 


25.63±0.45 


4.530 


12.25851 


03 32 48.89 


-27 52 43.17 


27.75±0.18 


26.49±0.12 


26.39±0.14 


>26.0 


>25.3 


>25.8 


>25.4 


4.540 


12.25952 


03 32 49.15 


-27 50 22.52 


28.24±0.29 


26.04±0.08 


25.30±0.05 


25.27±0.19 


25.52±0.39 


25.20±0.18 


25.54±0.33 


5.360 


12.26198 


03 32 49.81 


-27 50 22.75 


28.23±0.28 


26.42±0.10 


26.12±0.10 


25.58±0.27 


25.35±0.35 


25.94±0.41 


26.22±0.65 


5.110 


12.26409 


03 32 50.44 


-27 50 39.64 


>29.5 


27.01±0.20 


26.38±0.14 


>26.0 


>25.3 


26.30±0.57 


26.70±0.99 


5.460 


13.26480 


03 32 50.63 


-27 49 34.79 


>29.5 


26.83±0.16 


26.36±0.14 


>26.0 


>25.3 


>25.8 


>25.4 


5.220 


13.26492 


03 32 50.65 


-27 47 15.18 


26.80±0.09 


25.84±0.08 


25.64±0.07 


26.26±0.49 


>25.3 


25.42±0.28 


>25.4 


4.430 


12.26985 


03 32 51.94 


-27 52 8.494 


27.46±0.15 


26.04±0.09 


25.95±0.09 


26.19±0.58 


>25.3 


>25.8 


>25.4 


4.580 


12.27749 


03 32 54.05 


-27 51 12.02 


29.12±0.68 


27.09±0.20 


26.32±0.13 


>26.0 


>25.3 


25.82±0.85 


>25.4 


5.500 


12.28370 


03 32 56.22 


-27 51 51.29 


27.81±0.20 


26.33±0.10 


26.17±0.11 


>26.0 


>25.3 


25.01±0.21 


>25.4 


4.570 


12.28389 


03 32 56.29 


-27 53 31.53 


27.33±0.12 


25.66±0.05 


25.60±0.07 


24.62±0.17 


24.83±0.27 


24.61±0.15 


25.03±0.29 


4.690 


12.28728 


03 32 57.68 


-27 53 19.67 


28.05±0.25 


26.09±0.08 


26.08±0.10 


25.99±0.59 


>25.3 


24.99±0.21 


25.26±0.32 


4.930 


12.28859 


03 32 58.38 


-27 53 39.59 


26.41±0.05 


25.54±0.04 


25.62±0.06 


>26.0 


26.65±0.88 


24.98±0.30 


24.92±0.28 


4.420 


12.28917 


03 32 58.66 


-27 52 43.69 


28.33±0.32 


26.11±0.09 


25.81±0.08 


>26.0 


24.98±0.33 


25.05±0.20 


25.94±0.64 


4.840 


12.28990 


03 32 59.01 


-27 53 32.22 


27.47±0.13 


25.55±0.05 


25.20±0.04 


25.13±0.21 


24.42±0.17 


23.39±0.06 


24.07±0.13 


4.860 



to 
00 



Table 1 — Continued 



ID RA (,72000) Doc (J2000) v i' z J K^ 3.6 pm 4.5 /^m Zp^ot 

CO 

12.29097 03 32 59.72 -27 52 2.582 29.20±0.73 26.31±0.10 25.67±0.07 >26.0 >25.3 24.05±0.11 24.66±0.25 5.170 

12.29119 03 32 59.89 -27 52 56.42 28.53±0.40 26.70±0.15 26.29±0.13 >26.0 24.91±0.15 23.88±0.08 24.12±0.14 4.890 
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Table 2. Spectroscopically-confirmed z ~5 galaxies in GOODS-S Field 



ID 



RA (J2000) Dec (.12000) Rodshift v - i' 



Spitzcr Confusion Class 



Rodshift Flag 



Included in SED fitting 



44.1543 


03 32 5.258 


-27 43 0.406 


4.81 


1.85 


-0.09 


35.4142 


03 32 11.71 


-27 41 49.59 


4.91 


1.79 


0.17 


35.4244 


03 32 11.92 


-27 41 57.09 


5.57 


1.39 


1.05 


35.6626 


03 32 16.55 


-27 41 3.203 


5.25 


2.07 


0.92 


35.6867 


03 32 17.00 


-27 41 13.71 


4.41 


1.54 


0.11 


33.7471 


03 32 17.95 


-27 48 17.01 


5.40 


1.82 


1.16 


32.8020 


03 32 18.91 


-27 53 2.746 


5.55 


2.60 


0.57 


35.9350 


03 32 21.30 


-27 40 51.20 


5.29 


1.84 


0.71 


34.9738 


03 32 21.93 


-27 45 33.07 


4.78 


2.00 


0.28 


32.10232 


03 32 22.71 


-27 51 54.40 


4.90 


2.05 


0.42 


33.10340 


03 32 22.88 


-27 47 27.56 


4.44 


1.60 


0.00 


io.ACS.01 


03 32 22.89 


-27 45 20.99 


5.12 






33.10388 


03 32 22.97 


-27 46 29.09 


4.50 


1.65 


0.03 


34.11820 


03 32 25.31 


-27 45 30.85 


4.99 


3.55 


0.45 


35.14097 


03 32 28.56 


-27 40 55.71 


4.59 


1.67 


0.37 


35.14303 


03 32 28.84 


-27 41 32.70 


4.80 


1.79 


0.02 


io.ACS.02 


03 32 28.93 


-27 41 28.19 


4.88 






31.14602 


03 32 29.29 


-27 56 19.46 


4.76 


1.67 


0.10 


22.15184 


03 32 30.09 


-27 50 57.72 


5.08 


0.39 


-0.05 


24.18073 


03 32 34.48 


-27 44 3.008 


4.94 


1.46 


-0.12 


22.20159 


03 32 37.62 


-27 50 22.38 


5.51 


1.89 


1.28 


22.21502 


03 32 39.81 


-27 52 58.09 


5.54 


1.63 


1.03 


24.21686 


03 32 40.11 


-27 45 35.49 


4.77 


1.62 


0.16 


21.23040 


03 32 42.62 


-27 54 28.95 


4.40 


1.85 


0.41 


23.23051 


03 32 42.65 


-27 49 38.99 


4.84 


2.23 


0.06 


io.ACS.03 


03 32 43.15 


-27 50 34.80 


4.83 






23.24305 


03 32 45.23 


-27 49 9.829 


5.58 


1.61 


1.22 


21.24396 


03 32 45.42 


-27 54 38.52 


5.37 


2.64 


0.69 


22.25323 


03 32 47.58 


-27 52 28.18 


4.75 


2.01 


0.11 


12.28085 


03 32 55.08 


-27 54 14.48 


4.71 


1.88 


0.13 



4 
3 
3/4 
3 
1 
4 
1 
1 
3 
1 
1 

3/4 
2 
3 
4 



3/4 
3/4 

4 
3/4 

3 

3 

3/4 
3 
3 

3/4 



A 
C 
B 
C 
B 
C 
A 
A 
C 
C 
B 
C 
C 
B 
B 
B 
B 
B 
B 
C 
A 
C 
B 
C 
C 
C 
B 
A 
C 
A 



N 
Y 
N 
Y 
N 
N 
Y 
Y 
Y 
Y 
Y 
N 
N 
Y 
Y 
N 
N 
Y 
N 
N 
N 
N 
N 
Y 
Y 
N 
N 
Y 
Y 
N 



Note. — All magnitudes arc in AB system. no.ACS.01, no.ACS.02, and no.ACS.03, were not detected with ACS. The Spitzer 
confusion classes have the following meanings: 1 — isolated and detected; 2 — isolated but undetected; 3 — confused, but Galfit may 
help ; 4 — hopelessly confused. Those sources with Spitzer confusion class '3/4' were deemed hopelessly confused after attempting (and 
failing) to subtract nearby sources with Galfit. 
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Table 3. Photometric Properties of z~5 spectroscopically-confirmed galaxies. 



ID 


'^spec 


V 


i 


z' 


,7 


Ks 


3.6 Aim 


4.5 fini 


35.4142 


4.912 


27.22±0.11 


25.51±0.05 


25.26±0.05 


25.19±0.22 


25.09±0.35 


25.11±0.27 


25.65±0.54 


35.6626 


5.250 


29.07±0.62 


27.18±0.21 


26.35±0.13 






26.19±0.52 


>25.4 


35.6867 


4.416 


26.89±0.08 


25.38±0.04 


25.13±0.04 


24.92±0.23 


24.37±0.24 


23.43±0.04 


23.82±0.08 


32.8020 


5.554 


27.77±0.19 


25.13±0.03 


24.49±0.03 


24.74±0.13 


24.06±0.13 


22.73±0.02 


22.74±0.03 


35.9350 


5.283 


28.07±0.32 


26.04±0.10 


25.41±0.08 








25.52±0.60 


34.9738 


4.788 


28.22±0.29 


26.20±0.09 


25.82±0.09 


26.20±0.51 


24.90±0.23 


24.31±0.09 


24.72±0.18 


32.10232 


4.900 


27.90±0.25 


26.14±0.10 


25.68±0.08 


25.58±0.28 


25.03±0.24 


24.27±0.08 


24.82±0.16 


33.10340 


4.440 


26.64±0.07 


24.94±0.04 


24.84±0.04 


24.55±0.13 


24.59±0.16 


23.75±0.05 


24.01±0.10 


34.11820 


4.992 


28.78±0.45 


26.95±0.16 


26.66±0.16 


>26.0 


25.63±0.46 


>25.8 


>25.4 


35.14097 


4.597 


27.72±0.17 


25.92±0.07 


25.87±0.08 






25.19±0.30 


25.16±0.38 


31.14602 


4.760 


26.76±0.09 


25.04±0.04 


24.88±0.04 






22.64±0.02 


22.55±0.02 


21.23040 


4.400 


28.17±0.28 


26.12±0.08 


25.79±0.08 


25.30±0.26 


25.49±0.43 


24.02±0.06 


24.57±0.14 


23.23051 


4.840 


28.31±0.31 


26.04±0.08 


25.86±0.09 


26.19±0.38 


25.87±0.50 


24.93±0.14 


25.12±0.23 


21.24396 


5.370 


28.94±0.55 


26.09±0.08 


25.30±0.05 


24.99±0.22 


25.11±0.33 


24.80±0.16 


24.21±0.11 


22.25323 


4.758 


28.49±0.34 


26.43±0.11 


26.24±0.11 


25.43±0.26 


>25.3 


25.43±0.29 


>25.4 



Note. — The VLT mosaic does not cover the entire GOODS field. Those sources that are located off' the edge of the VLT images 
denoted by an ellipsis. 
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Table 4. Modeling Results 



Object Star Formation History log Mstciiar( Mq) Age (Myr) E(B-V) x^ 



35_4142 

35_6626 

35_6867 

32_8020 

35_9350 

34_9738 

32_10232 

33_10340 

35_14097 

31_14602 

21^3040 

23^3051 

21^4396 

22^5323 



csf 


9.34 


161 


0.00 


2.11 


70 


9.32 


143 


0.00 


0.26 


100 


10.37 


360 


0.01 


1.77 


300 


11.16 


905 


0.00 


2.90 


csf 


9.33 


255 


0.00 


0.57 


100 


10.13 


360 


0.00 


1.37 


70 


10.06 


255 


0.01 


2.22 


100 


11.28 


18 


0.24 


1.82 


200 


9.93 


255 


0.05 


0.04 


300 


11.10 


1015 


0.00 


1.79 


arst 


8.43 


1 


0.53 


2.24 


100 


9.86 


286 


0.00 


0.47 


arst 


8.40 


9 


0.17 


0.21 


arst 


8.43 


3 


0.32 


0.61 



Note. — In the star-formation history column, 'csf corresponds to constant star forma- 
tion, while the numbers (e.g 70, 100) correspond to the exponential decay constant (in Myr) 
for an exponentially-declining star formation history. 
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Fig. 1. — [Vqoqw — i775w) ^s. (^775vy — ^850Lp) colors of z~5 Candidates in GOODS-S. We 
construct a sample of 214 objects with photometric redshifts between 4.4 < z < 5.6 from 
the GOODS-MUSIC catalog (solid black circles). After removing stellar contaminants, low-z 
interlopers, and objects blended in Spitzer images, 72 objects with ^^ssolp < 26.5 remain; 
these objects are marked with an additional circle. Although ma ny objects in the sample 
fall just outside of the t>-band dropout selection window used by iGiavalisco et al.l (l2004al ) 
to select 2; ~ 5 galaxi es (demarcated by dotted line), redshift tracks generated from star- 
burst templates from iBruzual fc CharlotI (120031 ) illustrate that their rest-frame UV colors 
are consistent with the z~5 interpretation. These tracks assume an age of 100 Myr and 
a constant star formation rate with E(B—V) =0.0,0. 1,0.2 (blue dashed line, dashed-dotted, 



dashed-triple-dotted, respectively). Redshifting an elliptical galaxy template (jColeman et al. 



1980l ) to z=0-4 (red solid line), we see that old galaxies at z~1.5 could contaminate our sam- 
nle. 
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Fig. 2. — {vqoqw — ijY^w) ^s. (^775iy ~^850Lp) colors of 30 galaxies with F0R S2/VLT spectro- 
scopic r edshifts of 4.4 < z < 5.6. The f-drop selection window described in iGiavalisco et al. 
( l2004al ) is overlaid with a solid line. Starburst redshift tracks are identical to those described 



in Figure [TJ Removing objects that either do not satisfy the magnitude limit of ^^gsoLP <26.5 

or do not have the rest-frame UV colors expected for a z 

objects. 
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Fig. 3. — Distribution of IRAC 3.6 /iin AB magnitudes for 72 photometrically- selected z~5 
candidates (open circles) and 14 spectroscopically-confirmed z~5 galaxies (closed circles). 
The spectroscopic sample contains a larger relative fraction of Spitzer bright objects. 
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Fig. 4.— Observed and best-fit model iBruzual fc CharlotI fl2003f ) SEDs of 15 
spectroscopically-confirmed 2; ~5 galaxies in GOODS-S. Best-fit model parameters are pre- 
sented in Table H Three objects (32_8020, 31_14602, and 33_10340) have inferred stellar 
masses above 10^^ M©, and an additional three objects have inferred stellar masses greater 
than 10^° M0. 
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Fig. 5. — Confidence intervals for inferred stellar mass versus age for two objects from the 
spectroscopically-confirmed z~5 objects in GOODS-S. The ellipses are different assumed star 
formation histories, ranging from an initial burst to continuous star formation via a range 
of exponentially-decaying star formation histories. Contours are 68% confidence (solid line) 
and 95% confidence (dashed line) corresponding to ^xleduced = 1) 4 (respectively), where 
^Xreduced=xleduced ' xleduced,min- The Vertical dashed-dotted line at the right of each plot 
denotes the age of the universe at the source's redshift. Solutions to the right of this line are 
ruled out. The typical 68% confidence uncertainties in the stellar mass are are 30-50%. 
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Fig. 6. — Comoving star formation rate density as a function of redshift, assuming no 
extinction. The star formation rate densities are derived from Giavalisco et al. (2004) at 
z = 5, Stanway (2005) and Bouwens et al. (2006) aX z = 6, Bouwens et al. (2005) at 
z = 7.5, and Bouwens et al. (2005) at z = 10. The solid circles represent the star formation 
rate densities achieved by integrating the derived luminosity function down to 0.1 L^^g. 
Integrating the luminosity function down to zero luminosity (open circles) assuming a faint- 
end slope of q; = —1.73 adds an additional factor of 2.3 to the star formation rate density. 
The evolution of the star formation rate density with redshift is well fit by a (1 + z 
parameterization over 5 < z < 10. 
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Fig. 7. — Comparison of stellar mass density at z ^5-10 derived from spectral energy dis- 
tributions of galaxies in GOODS-S with that derived from integrating the observed star 
formation rate density at z^ 5 — 10 integrated down to 0.1 \j* z=z (solid line) and zero lu- 
minosity (dotted line) assuming a faint-end slope of a = —1.73. The filled solid circle with 
the upward arrow corresponds to the stellar mass density derived from spectroscopically- 
confirmed galaxies in GOODS-S (1x10^ Mq) whereas the solid filled circle corresponds to 
the stellar mass density from the photometric sample (6x10^ Mq). We estimate that the 
stellar mass density may be as high as 1x10^ M© Mpc~^ (solid horizontal line) depending on 
the contribution from undetected sources. We also include previous estimates of the stellar 
mass density at 2; > 5 from Yan et al. (2006) (diamonds), Eyles et al. (2006) (triangle), and 
Mobasher et al. (2005) (cross). The Yan et al. (2006) and Eyles et al. (2006) symbols are 
offset slightly from 2; = 6 for clarity. The large stellar mass density at z ~5 inferred from 
this study suggests that a significant amount of star formation is hidden by dust or has yet 
to be located, perhaps lying at higher redshift or in intrinsically faint systems. 



